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Panoramic Ha and mid-infrared mapping of star formation in a 
z = 0.8 cluster 
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ABSTRACT 

We present the first wide-field Ha imaging survey around the distant cluster RXJ 17 16.4+6708 
at z = 0.81 with a narrow-band filter on MOIRCS/Subaru, which reveals the star formation 
activities down to a star formation rate (SFR) of ^lM Q /yr without extinction correction. 
Combining with a wide-field mid-infrared (MIR) imaging survey with AKARI satellite, we 
compare in detail the unobscured and obscured star formation activities in the cluster. We 
find that both Ha emitters and MIR galaxies avoid the cluster central region and their spatial 
distribution is quite similar. Most of the Ha emitters show blue colours, but we find some 
Ha emitters on the red sequence. The MIR galaxies tend to be systematically redder than 
the Ha emitters probably due to heavy dust extinction. Interestingly, the red Ha emitters 
and the red MIR galaxies (i.e. dusty red galaxies) are most commonly seen in the medium- 
density environment such as cluster outskirts, groups and filaments, where optical colours of 
galaxies change. We investigate the amount of hidden star formation by calculating a ratio, 
SFR(IR)/SFR(Ha), and find that An a exceeds ~ 3 in extreme cases for actively star-forming 
galaxies with SFR(IR)>20Af Q /yr. It is notable that most of such very dusty galaxies with 
Au a >3 are also located in the medium-density environment. These findings suggest that dusty 
star formation is triggered in the in-fall region of the cluster, implying a probable link between 
galaxy transition and dusty star formation. We finally calculate the cluster total SFR and find 
that the cluster total SFR based on Ha alone can be underestimated more than factor ^2 even 
after 1 mag extinction correction. We suggest that the mass-normalized cluster SFR rapidly 
declines since z ~ 1 following oc (1 + z) 6 , although the uncertainty is still large. 
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1 INTRODUCTION 

It is well known that properties of galaxies are strongly correlated 
with environment in the sense that galaxies in high-density en- 
vironment tend to be red and have early-type morphology. This 
is first quantitatively noted by iDresslerl dl98d) in nearby clusters, 
and many other studies have confirmed such correlations as well 
as extended them in redshift space ( e.g. IPostman & Gelled fl98l: 
iDressler et al.ll997l : lGoto et al.l2003l ; |Postman et alj2005h . Recent 
large, intensive spectroscopic surveys for the local Universe show 
that star-forming activi t y is also a strong function of environm ent 
(e.g. iLewis et alj|2002l : [Gomez et al.ll2003l : iTanaka et alj|2004h in 
the sense that star-forming activity is weaker in the higher-density 
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environment. Cluster galaxies are in fact mostly red early-types 
with little on-going star formation at least in the local Universe. 
However, if we look at distant clusters, we can find many blue star- 
forming galaxi es in cluster environmen t (the so-called Butcher- 
Oemler effect: iButcher & Oemledll984h . Therefore by observing 
distant clusters at various redshifts (i.e. clusters in the past Uni- 
verse at various epochs) we should be able to identify when and 
how the activities of cluster galaxies are changed from blue star- 
forming populations to red quiescent ones. 

The 'hidden' star forming galaxies in clusters are also key 
population to characterize the evolution of cluster galaxies. Recent 
mid-infrared (MIR) observations of distant clusters with space tele- 
scopes such as Spitzer and AKARI discove red many dusty galax- 
ies in z £ 0-5 clusters dGeach et al.ll2006l; iMarcillac et al.ll2007l: 
iBai et alj [20071: ISaintonge, Tran. & HolderJ 120081; iKovama et al.l 
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120081: iKrick et al] 12009 ). Although the number of cluste rs stud- 
ied in MIR is still small. lSaintonge. Tran. & Holdenl ( 120081) showed 
the trend that the dusty star-forming populati on is more commonl y 
seen in the more di stant clusters (see also lHaines et al.l [2009b). 
iKovama et alj J2008t) showed the possibility of environmental de- 
pendence in the fraction of dusty star-forming galaxies at z ~ 0.8, 
and suggested that the MIR activity is enhanced in the 'medium- 
density' environment such as cluster outskirts, groups and fila- 
ments, where galaxy colours sharply change from blue to red. 

Another good tracer of star formation is optical emission 
lines from galaxies. In particular, Ha line (6563A) is considered 
to be one of the best star formation indicators, since it directly 
reflects the emission from Hll region in the star-forming sites. 
Moreover, Ha line is less affected by dust extinction or metal- 
licity, compared to [Oil] line which has been commonly used in 
the studies of distant galaxies. However, since Ha line shifts to 
near-infrared (NIR) regime at z > 0.4, it has been difficult to 
conduct a large Ha study for z <; 0.4 clusters until recently. 
There are some Ha studies for low-z clusters with optical spec- 
troscopy dCouch et alj 2001; Bal ogh et alj|2002[) and with narrow- 
band imaging dBalogh & Morrij2 000 ; Kod ama et alj2004l) , but the 
number of known clusters studied in Ha is r eally small. For distan t 
clusters. iFrnn. Zaritskv. & McCarthy! d2004l) and lFinn et all d2005h 
conducted narrow-band Ha imaging surveys for several z ~0.7- 
0.8 clusters in NIR. However, their field coverag e is limited to 
the ce ntral regions of the clusters. As suggested in lKovama et al.l 
(2008), cluster surrounding environment is likely to be the key en- 
vironment for galaxy evolution, but such spatially extended regions 
have been much less explored mainly due to the limited field of 
view of NIR instruments. 

It is thus ideal to observe distant clusters in both of the two 
major, but independent indicators of star formation, namely, Ha 
and MIR, over a wide area from cluster cores to surrounding envi- 
ronments. It is also important to quantify any different views of star 
forming activities seen by the two different indicators. Mapping out 
dusty and non-dusty star formation activities in distant clusters may 
give us a clue to understand what is actually happening in the tran- 
sition environments. In this paper, we present the first such wide- 
field Ha imaging survey over the known structures including the 
central cluster RXJ1716.4+6708 at z ~ 0.8. Combining with the 
similarly wide-field MIR data of this cluster taken with AKARI, 
we also attempt to compare the Ha and MIR views of the z > 0.8 
cluster. 

The RXJ1716.4+6708 cluster (hereafter RXJ1716) that we 
target in this paper is a rich and probably unvirialized cluster 
at z — 0.81. The cluster was first discovered i n the ROSAT 
North Ecliptic Pole Survey (NEP: Henrv et al. 1997). Optical spec- 
troscopy was performed bv lGioia et alj d 19991) and 37 c luster mem- 
bers w ere identified. Using this spectroscopic sample, iGioia et all 
( 1 19991) determined the cluster redshift z c i = 0.809 and the ve- 
locity dispersion a = 1522+jgQ km s" 1 . This velocity dispersion 
is relatively large for its rest-frame X-ray luminosity of Ltoi = 
13.86±1.04xl0 4 4 erg s" 1 and the temperature kT = 6.8^ keV 
dEttori et alj|2004l) . This indicates that this cluster is not virialized 
yet and in the process of active assembly. In fact, this cluster has a 
small subcluster or group to the northeast of the main cluster. The 
morphology of this cluster in a X-ray image elonga t es towards the 
directi on of the subcluster (e.g. ljeltema et alJl2005h . lKovama et al.1 
d2007I) performed a wide-field, multi-colour optical imaging of this 
cluster and discovered prominent large-scale structures penetrating 
the cluster core and the second group of this cluster towards the 
southwest of the cluster core, based on the photometric redshift 



technique. The weak-len sing mass of this c luster is estimated to be 
2.6±0.9xl0 14 /i~ 1 Af Q dClowe et al.ll998l) . T his is consist e nt wit h 
the estimated mass based on the X-ray data in lEttori et al.1 i2004) : 
Mtot = 4.35 ± 0.83 x 10 14 M Q . 

The structure of this paper is as follows. In Section 2, we 
summarize our optical, NIR and MIR data. In Section 3, we show 
the selection technique of Ha emitters from cluster members. The 
derivation of Ha-derived star-formation rates is shown in Section 4. 
We show the results and discussions from Section 5 to Section 8, 
and summarize our results in Section 9. Throughout this paper, we 
use fW = 0.3, Oa = 0.7, and H = 70 km s _1 Mpc -1 . Magni- 
tudes are all given in the AB system. 



2 DATA 

2.1 Optical data 

We use the optical data in IKovama et al.1 d2007l) . We performed 
deep and wide-field VRi 1 V imaging of the RXJ1716 cluster 
with Suprime-Cam dMivazaki et al.1 |2002) on the Subaru Tele- 
scope dive et al.ll2004h . as a part of the PISCES project (Panoramic 
Imaging and Spect roscopy of Cluster Evolution with Subaru; 
Koda ma et al .1120051) . A summ ary of the data and t he method of 
data reduction is described in IKovama et alj d2007l) . and we just 
repeat some important points here. 

We reduced the data using a pipeline software SDFRED 
dYagi et al.1 12002| ; lOuchi et al.1 |2004|) . Source d etection and pho- 
tomet ry were done using SExtractor software dBertin & Amoutsl 
1996). We use 2" aperture magnitudes (MAG.APER) for mea- 
suring galaxy colours and MAG_AUTO are used as total magni- 
tudes of galaxies. Our optical catalogue was constructed for ob- 
jects brighter than z = 24.9 which corresponds to the 5a detec- 
tion limit in z'-band. For our use in this paper, we applied a small 
aperture correction to the aperture magnitudes. This is because the 
seeing size of our near-infrared (NIR) data is slightly larger (0.89", 
see Section 2.2.1) than that of the optical images (~ 0.7"). We 
smoothed our optical images to 0.89" and determined the correc- 
tion value, but the correction is small (~0.02 mag). We estimated 
photom etric redshifts (phot-z) fo r all th e galaxies using the phot-z 
code of lKodama. Bell. & Bowen dl999h . We fo und that > 90% of 
spectroscopically confirmed members listed in IGioia et al.1 ( fl999h 
fall within the red shift range of 0.76 < z P hoi < 0.83 (see fig. 2 of 
IKovama et alj|2007l) , and used this range to select member candi- 
dates. 



2.2 Near-infrared(NIR) data 

2.2.1 Observation and data reduction 

We observed the RXJ1716 cluster in J and NB119 filter 
with MOIRCS on the Subaru Telescope dlchikawa etal] 120061; 
ISuzuki et alj2008h . The NB 1 19 filter perfectly matches to Ha lines 
at z = 0.81 (i.e., the redshift of our target cluster). In Fig. 1 we 
compare the veloc ity distribution of spectroscopically confirmed 
members listed in IGioia et al", I d 19991) with the NB119 filter re- 
sponse function. It clearly shows that we can detect more than 
~80% of Ha emitting members within the FWHM of the NB119 
filter if they have the s ame velocity distr ibution as that of the spec- 
troscopic members in lGioia et lji1( fl999T) and if their Ha EWs are 
larger than the threshold of our observation. The J-band filter sam- 
ples mainly continuum flux. Therefore, with a combination of J 
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Figure 1. The transmission curve of the NB1 19 filter. The histogram shows 
the radial velocity distribution of spectroscopically confirmed members 
listed in Gioia et al. (1999), with respect to the velocity centre of the clus- 
ter. The wavelength of the histogram indicates the position where Ha line 
of each member galaxy falls. 



and NB119 filters, we can conduct an unbiased, deep Ha imag- 
ing survey of this cluster. To neatl y cover the large-scale structures 
discovered bv lKovamae7ai] (2007), we set eight fields of view of 
MOIRCS as shown in Fig. 2. It can be seen that Fl, F2, F7 and F8 
have the full instrument coverage of 4' x7' (i.e. the original FoV of 
MOIRCS), and all the other fields, F3, F4, F5 and F6, have only 
a half size. Unfortunately, at the time of observation, one of the 
two MOIRCS chips was replaced with an engineering-grade one 
due to the instrumental failure. We should still stress, however, that 
this is the first Ha imaging survey covering such a wide range in 
environment at z > 0.8. 

The data is reduced in a standard manner primarily using a 
pipeline software MCSRED (Tanaka et al. in prep). Some frames 
taken with NB1 19 show a strong fringe pattern, and we used a self- 
made fringe subtraction software. After matching the source posi- 
tions on individual frames, we coadded to make a final image for 
each field, and then mosaiced all the images together. A summary 
of the J and the NB119 data is shown in Table 1. We match the 
seeing size of all the images to 0.89", which is the worst seeing 
size in the J-band. The photometric zero-points are derived using 
some standard stars observed on 27th May 2007. We properly scale 
F3-F8 data to Fl or F2 using the photometry of the same objects 
in the overlapping regions, although the amount of such scaling 
is small (<0.1 mag). We finally check that there is no significant 
difference in the z ero-points between our FoVs using stars in the 
2MASS catalogue Jjarrett et ai1l2000l) that fall within our observed 
fields. The Galactic extinction is cor rected using the dust map by 
Schlegel, Finkbe iner, & Davisl (l998). Limiting magnitudes are es- 
timated by measuring the scatter in fluxes in randomly distributed 
2" apertures for each FoV. Note that some overlapped regions are 
deeper. 



2.2.2 Photometric catalogue 

We first co nstruct an object cata logue in J and NB119 with the 
SExtractor ( Berti n & Arnoutsl 19961) using its two-image mode with 
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Figure 2. Our MOIRCS FoVs (F1-F8) are shown by the solid-line boxes. 
The three large dashed-line squares indicate the FoVs of our MIR observa- 
tions with AKARI satellite (10' X 10' for each). The phot- ,2 selected mem- 
ber galaxy candidates in Koyama et al. (2007) are shown by the small dots. 
The contours indicate l,2,3,4,5cr above the average local number density of 
the phot-z selected member candidates. 



the NB119 image for source detection. We use MAG_APER (2" 
aperture magnitude) for colours of galaxies and MAG_AUTO for 
total magnitudes. As shown in Table. 1, the depths of the data vary 
from field to field. We include objects with NB119 < 23.2, which 
is > 5cr detection in all the fields. We also require J < 23.9 
mag (~3<j in J-band) for our sample to avoid spurious detection 
in NB119. Af t er ma tching to the optical catalogue constructed by 
IKovama et al. (2007), we construct a J- and NB 1 19-detected cata- 
logue of the RXI1716 cluster. Stars and objects around very bright 
stars are rejected through this process. The catalogue contains a 
full range of information of the NB 1 19-detected galaxies, including 
photometrie s in VRi'z'J and N B119 bands as well as the phot-z 
values from Koya ma et al.l ( 120071) . 



2.3 Mid-infrared(MIR) data 

iKovama et alj d2008l) observed this cluster at 15/im (LI 5 fil- 
ter) with infrared camera (IRC) onboar d the AKARI satellite 
dOnaka et alj|200l iMurakami et al.ll2007l) . The FoVs of our MIR 
observations are shown with dashed-line squares in Fig. [2] It is 
notable that the areas of our NIR observations are mostly cov- 
ered by those of MIR observations. Our 15/im fil ter captures the 
peaks of Polycyclic Aromatic Hydrocarbon (PAH: IPuget & Legej 
Il983) emissions at rest-frame 7.7 and 8.6 /iin at z ~ 0.8, which 
are considered to be a good indicator of dusty star formation (e.g. 
ICharv&Elbazll200lh . The observed 15 /im flux (FLUX_AUTO 
value from SExtractor) is first converted to the total IR lumi- 
nosity, i(IR), using a relation b etween vL v ,g, um and L(IR) de- 
rived from the starburst SEDs in lLagache et ail (2004) (see fig.4 
of IKovama etai]|2008l) . Thus derived L(IR) is then converted to 
SFR using the c onversion factor between SFR and L(IR) from 
Kennicutt (1998), i.e., SFR[A/ Q /yr] = 4.5 x 10~ 44 L(IR)[erg/s]. 
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Table 1. A summary of our near-infrared imaging data in J and NB 119. 



Note that our 5<r detection limit is 67^iJy which corresponds to 
SFR(IR)~15 Af /yr. We adopt a 50% uncertainty in SFR(IR), 
which is larger than simple photometric errors, to take into ac- 
count the rep orted uncertainty i n the conversion fro m vL^^^m 
to L(IR) (e.g. lCaputi et al.ll2007l: iBavouzet et aj]|2008f) . Using the 
MIR catalogue constructed by Ikovama et al.l d2008f) 7we made a 
catalogue of MIR-dete cted galaxies within th e fields covered by our 
MOIRCS observation. IKovama et af] d2008h constructed two types 
of catalogues; one for the 15/im-resolved members (single objects) 
and the other for the 15/im-unresolved members (blended objects). 
Likewise, we make such two types of MIR-detected catalogues in 
this paper as well. We derive star formation rates in IR, SFR(IR), 
for the 15/im-resolved members only, since there are large uncer- 
tainty in the MIR photometries of the blended sources. 




18 19 20 21 22 23 
NB119 [mag] 



3 SELECTION OF Ha EMITTERS FROM CLUSTER 
MEMBERS 

3.1 Selection technique 

If a line emission comes into the NB119 filter, it can be identi- 
fied as an object with a flux excess in NB119 compared to the J- 
band flux. Therefore the J— NB119 colour is a good indicator to 
search for emitters. However, we should take into account the fact 
that the NB119 filter (A c = 11885A) is located near the lower- 
wavelength end of J-band (A c = 12600A). Therefore, we cannot 
simply use the J-band magnitude as the continuum flux because 
the continuum flux at 1.19 pirn and the J-band flux can be differ- 
ent depending on the SED slope around that wavelength regime 
of each galaxy. Therefore, we estimate the continuum flux level at 
1.19/im of each galaxy by linearly interpolating the magnitudes in 
z and J bands. This correction substantially reduces the scatter 
in J— NB119 colours around the zero value which is expected for 
non-emitters in NB119 at various redshifts. We refer to such flux 
'corrected' J-band magnitude as J CO rr, and J CO rr — NB119 colours 
of our sample are plotted against NB 1 19 magnitudes in Fig. 3. The 
solid curves show ±3cr errors in J— NB119 colours. We define 
NB119 emitters as those having colour excesses of J corr — NB119 
> 0.3 and above the upper curve at the same time. 

The NB119 emitters thus selected are not necessarily Ha 
emitters at 2 = 0.81 because other line emitters at other redshifts 
are also included in the NB119 emitters. For example, H/3/[OIII] 
emitters at 2 — 1.4-1.5 and [Oil] emitters at 2 — 2.2 are the 
major contaminants. We eliminate these galaxies using broad-band 
colours. We show in Fig.|4]a colour-colour diagram. We plot phot- 
2 selected galaxies (i.e. 0.76< 2 p hot <0.83) as small filled circles 
and all the NB1 19 emitters as crosses. We also show the predicted 



Figure 3. The narrow-broad bands colour-magnitude diagram to define 
our NB119 emitters. The J-band magnitudes are corrected for a colour- 
term using z'~J colours so that they correspond to the continuum fluxes at 
1.19/im. The vertical dashed-line shows 5<r limiting magnitudes in NB119 
and the solid curves indicate ±3cr excess in J— NB119 colours. Galaxies 
above the horizontal solid line (J— NB119 =0.3) and the upper solid curve 
are denned as the NB119 emitters (shown by crosses), of which the open 
squares show the Ha emitters at z = 0.81 as defined later (see text and 
Fig.|4). 



colour tracks of model galaxies at 2 = 0.8, 1.45, 2.2, respectively, 
along which the contributio n of bulge component is chan ged from 
0.0 to 1.0 at each redshift jKodama, Bell, & Bowerll 19991) . We can 
see that some emitters have very different colours from the model 
track of 2 = 0.8 galaxies or the distribution of the phot-2 selected 
galaxies. These galaxies are likely to be contaminants. We set a 
boundary shown by a closed box in Fig. [4] and define the Ha emit- 
ters at 2 = 0.81 as the NB119 emitters that satisfy these criteria. 
We may miss a few real members near the boundary, but their effect 
to our conclusion is negligible. 

We summarize our definitions of cluster members in the fol- 
lowing. Firstly, we select all galaxies with 0.76 < z p hot < 0.83 
as members. This criterion recovers mor e than 90 % of sp ectro- 
scopical ly confirmed m embers listed in iGioia et alj d 19991) (see 
fi g. 2 i n Koyama et al. 2007 ). The same criterion is actually used 



in IKovama et alj d2007l) and IKovama et alj |2008). The number of 
these phot-2 selected members is 447. Secondly, we include all the 
Ha emitters at 2 = 0.81 selected above. We have 114 Ha emit- 
ters in total. Note that 32 of them have been missed out by our 
phot-2 selection, although many of them have 2 p hot ~ 0.70-0.75, 
only slightly lower than our phot-2 selection criterion. This is in- 



© 0000 RAS, MNRAS 000, 000-000 



A Ha and MIR imaging of a z = 0.8 cluster 5 



o 

-0.5 





1 1 1 1 | 1 1 1 1 | 1 1 1 1 | 1 1 
E(B-V)=0.2 


1 1 1 _ 

z=1.45 : 

H / 


- 


• J' ^'Jr>iv' 
1 = 0.8/ '.jSj^/ 

m » / 


»» / / - 

/ M / 
/ 8 / H 

/ " / z=2.2 : 

s « 










" »" 




- 





" 



0.5 1 1.5 2 2.5 3 
i'-J [mag] 



Figure 4. The colour-colour diagram to select Ha emitters at z = 0.81. 
The phot-z selected members (filled circles) and the NB119 emitters 
(crosses) are plotted. Three solid curves show the colour tracks of galaxies 
at 2 =0.8, 1.45, 2.2, respectively, based on the Kodama et al. (1999) model. 
Along the curves, the bulge fraction to the total light is changed from 0.0 to 
1.0 (blue side to red side) at fixed redshift. The NB119 emitters located in 
the closed box are defined as Ha emitters at z = 0.81. The arrow shows a 
reddening vector of E(B — V)=0.2, calculated from the extinction law of 
Calzetti et al. (2000) assuming R v =3.l. 
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Figure 5. Ha luminosity function of the RXJ1716 cluster (filled diamonds 
with purely poissonian error bars). The normalization of the cluster LF is 
arbitrarily scaled to compare with the field Ha luminosity function of So- 
bral et al. (2009) (filled circles and their best-fitted LF as shown by the solid 
curve). 

We note that as shown in iKovama et alj J2008h the optical colour 
distribution starts to change dramatically from blu e to red in the 
medium-density environm ent thus defined (see also lKodama et all 
l200ll:lTanaka et alj|2005h . 



evitable since the phot-z tends to give lower redshift estimation 
for star-forming galaxies due to the lack of prominent features in 
their SED, while redshift estimatio n for passively evolving galaxie s 
are relatively more accurate (e.g. iKodama. Bell. & Bower] 1 19991) . 
This means that the combined approach that we take in this paper, 
namely the phot-z selection and the narrow-band emitter survey, is 
quite effective to increase the completeness of true cluster mem- 
bers. 



3.2 Local density measurements 

One of the main purposes of this paper is to examine environmen- 
tal dependence of Ha-derived star formation activity at z = 0.8. 
To quantify the environment, we estimate the local density of each 
member galaxy using its five nearest neigh bours. We already de - 
fined the local density (£5 in Mpc -2 unit) in lKovama et al.l (2008). 
but we recalculated it since we added 32 new members in this pa- 
per as shown in the previous sectio n. When we ca l culate £5, we 
use all the member galaxies used in lKovama et al.1 J2OO8) (i.e. all 
the ~ 2700 galaxies throughout the Suprime-Cam field in the op- 
tical catalogue with 0.76 < z ph ot < 0.83 to the depth of z'=24.9) 
and the 32 newly selected members. Theref ore, it makes almost n o 
change for the estimation of £5 from that o f lKovama et al] (2008), 
but the main purpose here is to define environment with £5 in the 

same way for the newly added me mbers. 

Following the definition in Kov ama et al.l |2008), we clas- 
sify galaxies into three environmental bins, namely, 'low-density', 
'medium-density' and 'high-density' regions which correspond to 
log £5 < 1.65, 1.65 < log £5 < 2.15 and log £ 5 > 2.15, respec- 
tively. The high-density region corresponds to the cluster core. The 
medium-density region corresponds to cluster outskirts, groups, 
and filaments. The low-density region corresponds to outer fields. 



4 Ha-DERIVED STAR-FORMATION RATE AND 
SPECIFIC STAR FORMATION RATE 

4.1 Derivation of SFR(Ha) and SSFR(Ha) 

In this subsection, we calculate some basic physical quantities to 
characterize the star formation activity, such as star-formation rates 
(SFRs) and specific star-formation rates (SSFRs), for all the Ha 
emitters. Firstly, we derive Ha line flux, continuum flux and rest- 
frame equivalent width by the following equations: 



Ha + [NII] 



/nb - fj 
' 1 - Anb/Aj 



fj -/nb(Anb/Aj) 



1 - Anb/Aj 



EW fl (Ha+ [Nil]) = (l + z) 



-1 Fua 



(i) 



(2) 



(3) 



where Aj and Anb are the widths in wavelength of the J and 
NB119 filters, fj is the flux density at 1.19 /im derived from the 
J COT r, and /nb is the flux density for NB119. We then multiply 
47rd| to -Fhq+inii] to derive the luminosity L(Ha+ [Nil]), where 
d,L is the luminosity distance 5.10 x 10 3 Mpc at z = 0.81. Fi- 
nally, we compute t he Hq-based sta r formation rates, SFR(Ha), 
using the relation of lKennicuttl dl998l) , SFR(Ha)[M /yr]= 7.9 x 
10~ 42 LHa,[erg/s]. We note that our selection criteria of Ha emit- 
ters shown in Section 3.1 correspond to EWR(Ha+[NII]) >30A and 
SFR(Ha)^lA/0/yr before extin ction correction. W e correct for 
30% [Nil] line contributi on (e.g. Tresse et alJI 19991) and conside r 
1 mag extinction of Ha (Kennicutt, Tambl vn. & Congdonlll994b . 
We should note that these assumptions adopted here are somewhat 
uncertain and can only be applied to discuss statistical behaviours. 
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Figure 6. The physical quantities of the Ha emitters, EW(Ha+[NII]), SFR(Ha) and SSFR(Ha), are plotted as a function of the local density (left) and stellar 
mass (right). The open circles show all the Ha emitters, the filled circles represent the MIR galaxies and the open pentagons indicate the red Ha emitters 
with R — J > 2.0. In SFR and SSFR, we have applied a correction for 1 mag extinction. In the top panels, the non-Ha emitters are also plotted at EW=0, 
while only the Ha emitters are plotted in the middle and the lower panels. The two small arrows with numbers in the top panels indicate the galaxies with high 
EW(Ha+[NII]) located above the top boundary. In the left panels, we show median values as crosses in each environment with error bars showing 25-75% 
distribution of each sample. The vertical dashed lines at log £5 =1.65 and 2.15 in the left panels indicate the dividing lines between the low-, medium- and 
high-density environments as defined in Koyama et al. (2008). 



In particular, the amount of extinct ion can be much larg er for very 
actively star-forming galaxies (e.g. lPoggianti & WuluOOOl) and this 
will be discussed in Section 7 of this paper. We also note that we do 
not correct for the contamination of active galactic nuclei (AGNs) 
in our Ha emitters sample. In reality, some Ha emission lines 
may be originated from AGNs rather than star-forming Hll regions. 
However, it is reported that the fraction of AGNs in distant clusters 
is very small (only ~ 1%) from optical spectroscopic surveys (e.g. 
iDressler et a0l 999). The fraction of X-ray selected AGNs could be 
higher by a factor o f ~ 5 than that of optically selected ones (e.g. 
Martin i et ai]|2002l) . But in any case, AGN contamination is not a 
major concern and has li ttle effect on our statistical discussion on 
the Ha emitters (see also lGarn et al.|[2009h . 

One of the best indicators to quantify the status of galaxy ac- 
tivity is the specific star-formation rate (SSFR), which is SFR di- 
vided by stellar mass of galaxies (i.e. SFR per unit stellar mass). We 
derive stellar mass (A/») of galaxies in the following way. We use 
the relation between A/„/Lj( obs ) of galaxies and R — J at z — 0.8 
for the mod el galaxies with different contri bution of bulge com- 
ponent from lKodama. Bell. & Bowed (1999). Note that the stellar 
mass in the model is re-scaled using the Salpeter IMF for consis- 
tency with the derivation of SFR. Based on the observed R — J 
colour of each galaxy, we derive M*/Lj( obs ) which is then multi- 
plied to the observed J-band luminosity to get M* . SSFR is derived 
by dividing SFR by M*. 



4.2 Ha luminosity function 

In Fig. [5] we show the Ha luminosity function (LF) using Ha emit- 
ters within 2.5 Mpc from the cluster centre. We here slightly modify 
the selection criteria of Ha emitters to make a fair comparison with 



ISobral et al.l d2009h . who constructed a field Ha LF at z ~ 0.8. 
Assuming that all Ha emitters are associated with the cluster, we 
do not apply field subtractio n. Our cluster LF is arbitrarily scaled to 
match the field Ha LF from lSobral et al.l d2009l) (filled circles with 
their best-fitted LF). Therefore, we can only compare the shape of 
the cluster and field Ha LFs. It is apparent that there is no signif- 
icant differ ence in the shape of the Ha LF. This is similar to the 
result from Kodama et al. ( 2004). They showed that the Ha lumi- 
nosity function shows little sensitivity to local density, using their 
very wide Ha imaging around C L0024 cluster at z = 0.39. In the 
local U niverse, this is also true jBalogh et alj|2004h . iBalogh et al.l 
d2004> limit their sample for star-forming galaxies and showed that 
the strength of Ha lines does not depend on environment. Thus, al- 
though the fraction of star-forming galaxies may change as a func- 
tion of environment, star-forming activity of star-forming galaxies 
seems not strongly dependent on environment. Our new Ha study 
for the z — 0.81 cluster extends this idea up to z ~ 0.8. As will be 
discussed in Section 7, there are some exceptionally active galax- 
ies in the cluster surrounding environment. However, it would be 
difficult to find such excess of activity in this Ha LF. At least, our 
comparison suggests that we will need much larger sample to de- 
tect the environmental difference of star-forming activity among 
the star-forming galaxies, if any. 



4.3 Dependence on environment and stellar mass 

InFig.Hl we plot EW R (Ha+[NII]), SFR(Ha) and SSFR(Ha) of Ha 
emitters as a function of local density and stellar mass. It is apparent 
that very strong Ha emitters with EW>100A are seen only in the 
medium- and the low-density environments, and no such galaxy is 
found in the high-density environment (see the top-left panel of Fig. 



© 0000 RAS, MNRAS 000, 000-000 



A Ha and MIR imaging of a z = 0.8 cluster 7 



6), although this may be due to small number statistics of the Ha 
emitters in the high-density region. Other than that, we see no sig- 
nificant environmental trends in these quantities of the Ha emitters. 
Rather, they show modest correlations with stellar mass of galax- 
ies (right panels). Strong Ha emitters with EW^IOOA are found 
only among low mass galaxies with log M* < 10.4. More mas- 
sive galaxies have higher SFRs but lower SSFRs than less massive 
galaxies. 

We show in Fig. [7] the distribution of SSFR(Ha) for each 
environment, dividing our Ha emitter samples into low-mass 
(log M„/M©< 10.4) and high-mass (log M»/M Q > 10.4) galax- 
ies. We find that high-mass galaxies tend to have low SSFR 
(SSFR< — 9.5) compared to low-mass galaxies in all environment. 
The lack of low SSFR galaxies in low-mass galaxies might be 
partly due to the selection effect, since the selection criteria of NB 
emitters is more strict for fainter galaxies as can be seen in Fig. [3] 
However, the lack of high mass galaxies with high SSFR should be 
reliable. We show the median value of SSFR(Ha) in each panel as 
an arrow. We find that the median values do not strongly correlated 
with environment if we fix the stellar mass of galaxies, although 
the number of low-mass Ha emitters in high-density environment 
is really small. We should note that this suggestion is valid only 
when we limit the sample to the Ha emitters, and we should keep 
in mind that there are many Ha-undetected galaxies in particular in 
higher-density environments (see Section 5). 

We plot the MIR-detected galaxies as filled symbols in Fig. [6] 
We can notice that the MIR-detected galaxies tend to have low 
EW(Ha+[NII]) and low SSFR(Ha) in spite of their high star for- 
mation rates in MIR (SFR(IR);> 15M /yr). One may claim that 
MIR-detected galaxies tend to be biased to more massive galaxies 
with larger SFR. However, interestingly, there is no significant dif- 
ference in SFR(Ha) between MIR-detected and MIR-undetected 
galaxies (see middle-right panel in Fig. [6}. There is a trend that 
massive Ha emitters tend to be detected in MIR, indicating that 
they tend to be more dusty, hence their SFRs (and SSFRs) are likely 
to be underestimated even with Ha. This suggests that we need not 
only Ha information but also MIR information to fully understand 
the star formation history of massive galaxies in the distant Uni- 
verse. In the remaining of this paper, taking the great advantage of 
wide-field coverage of our Ha and MIR data, we attempt to show 
how large amount of star-formation is obscured and hidden in the 
previous optical surveys and how the obscuration is related to envi- 
ronment at 2 ~ 0.8 for the first time. 



5 MAPPING OUT STAR-FORMATION ACTIVITY IN 
AND AROUND THE CLUSTER 

5.1 Spatial distribution 

We show in the left panel of Fig. [8] the spatial distribution of clus- 
ter member galaxies (small open circles), blue Ha emitters with 
R — J <2.0 (open squares) and red Ha emitters with R — J >2.0 
(open pentagons). We only plot the galaxies in the area observed by 
the NB filter. We can see that Ha emitters well trace the filamen- 
tary large-scale structures. This is a strong evidence for that most 
of the small groups identified by phot-z are physically associated 
with the cluster and not just a concentration of fore-/background 
galaxies. Note that we may be missing some Ha emitters especially 
in the outskirt of the clusters, where relative velocities of galaxies 
with respect to the cluster core would be larger, although the width 
of the NB119 filter neatly covers the radial velocity range of ± 
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Figure 7. The distributions of SSFR for high-mass (log M*/M Q > 10.4) 
and low-mass (log M* /Mq < 10.4) galaxies in each environmental bin for 
all the Ha emitters. The arrows indicate the median value of SSFR in each 
panel. The vertical dashed lines at SSFR=— 9.5 show the dividing line be- 
tween the high and low SSFR galaxies for an eye-guide. 



~1500 km s _1 (see Fig. 1). A 5'x5' close-up view is shown in 
the right panel of Fig. [8] We also plot the position of MIR-detected 
galaxies (including resolved and unresolved objects) as filled cir- 
cles. It can be easily noticed that the distribution of Ha emitters 
is quite similar to that of MIR-detected galaxies. We reported the 
av oidance of MIR-dete cted galaxies in the cluster central region 
in lKovama et al.l ( 120080 , and this study strengthens the conclusion 
that the star-forming activity is very low in the cluster core of this 
cluster. 

On the other hand, it is interesting to note that the Ha emit- 
ters and MIR-detected galaxies do not always directly overlap each 
other. Since our Ha survey can go deeper in terms of SFR (~ 1 
M©/yr without extinction correction) compared to MIR observa- 
tion (~ 15 Mq/jt), it is natural to expect that some moderately 
(or weakly) star-forming galaxies are detected only in Ha. In fact, 
we can see many such galaxies in Fig. [8] However, there are also 
some galaxies which are detected in MIR but not in Ha. These type 
of galaxies are difficult to interpret. It may be possible that these 
galaxies are very dusty and a large amount of UV/optical light is 
hidden by dust, so that we cannot detect their Ha emission lines. 
It is also possible that some MIR galaxies are located outside of 
the redshift range of the Ha emitters that is covered by the NB1 19 
filter. Some of them can still be members (see Fig. 1), but some 
could be non-member contaminants due to the limitation of phot- 
z based membership determination of the MIR galaxies. Spectro- 
scopic confirmation of their membership is needed to discuss this 
interesting population further. 



5.2 Ha fraction 

As shown in the previous section, Ha emitters avoid cluster central 
region and/or high-density region. We quantify it by calculating 
the fraction of Ha emitter as functions of local density and cluster 
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Figure 8. (left): The 2-D distributions of all the members (open circles), the MIR-detected galaxies (filled circles), the blue Ha emitters with R — J<2.0 
(open squares), and the red Ha emitters with R — J>2.0 (open pentagons). The dashed-line boxes indicate our MOIRCS FoVs. We plot the member galaxies 
detected in our MOIRCS FoVs only. The contours are the same as in Fig. 2. (right): A close-up view of the cluster central region. The meanings of the 
symbols are the same as in the left panel, except that the large and the small open symbols indicate the strong (EW> 5()A) and weak (EW< 50A) Ha emitters, 
respectively. The solid-line circle shows the radius of 0.5 X R200 from the cluster centre (see Section 8). 



centric radius (Fig.[9]l. We use galaxies with NB119<22.0 where 
the selection of Ha emitters is complete down to J— NB 119=0.3 
(see Fig. [3}- It can be seen that the Ha fraction is strongly depen- 
dent on local density and cluster centric radius, in the sense that 
the fraction is lower in the higher-density environment. Thus, we 
conclude that there is little star-formation activity in the core of the 
RXJ1716 cluster. Such global trend that the fraction of Ha emit- 
ters decreases towards higher density regions (Fig. 9) is largely due 
to the fact that the fraction of passive galaxies increases towards 
higher density regions. In fact, if we limit the sample only to blue 
galaxies with R — J < 2.0, the fraction of Ha emitters does no 
longer strongly correlate with local density or cluster centric ra- 
dius, and is almost constant at ~60%. This is consistent with the 
results obtained in Section 4.2 and 4.3 that the observed Ha line 
strength in star-forming galaxies does not strongly correlate with 
environment. 

Finn ct al. I J2005h conducted similar Ha emitter survey for the 
central part of three z ~ 0.7 clusters. They showed that two of 
the three clusters have few Ha emitters in the central region. How- 
ever, th ey found many Ha e mitters in the core region of the other 
cluster. lHavashi et alj d2009h show that there are many [Oil] emit- 
ters even in the cluster core at z = 1.46, and that the fraction of 
[Oil] emitters is still high in the core within _R C <0.25 Mpc from the 
cluster centre. For the RXJ1716 cluster, we find no Ha emitters at 
7?c<0.25 Mpc as shown in the right-panel of Fig.|9] although there 
are 31 member galaxies in i? c <0.25 Mpc without significant detec- 
tion in Ha. Therefore, for the RXI1716 cluster, the termination of 
star-forming activity in the cluster core has been completed before 
z — 0.8. These facts suggest that the epoch of the termination of 
star-forming activity in the cluster core may differ from cluster to 
cluster, probably depending on the cluster mass and/or maturity of 
clusters. We need a much larger sample of distant clusters to draw 
a general picture of formation and evolution of the cluster cores. 




3 2.5 2 1.5 1 
log2 5 [Mpc-2] 



0.1 1 
Rc [Mpc] 



Figure 9. The fraction of the Ha emitters as a function of local density 
(left) and as a function of cluster-centric distance (right). The error-bars 
show Poisson errors. 



6 PROPERTIES OF Ha EMITTERS AND MIR 
GALAXIES 

6.1 Colour-magnitude diagram 

In Fig. [TO] we show (R — J) v.s. J colour-magnitude diagram 
using all the member galaxies. Ha emitters are marked with 
open squares again. The small and large sizes of squares indicate 
EW(Ha+[NII])<50A and EW(Ha+[NII])>50A, respectively. We 
clearly notice that Ha emitters are mainly blue galaxies and that 
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Figure 10. (left): The colour-magnitude diagram for all the observed fields. The open circles show all the member galaxies. The small and large open squares 
indicate the weak (EW<50A) and strong (EW>5()A) Ha emitters, respectively. The filled circles indicate the MIR galaxies. The histograms show the colour 
distributions of the Ha emitters and the MIR galaxies as indicated. The arrow shows a reddening vector that corresponds to E(B — V) = 0.2. (right): The 
same plot for each environment (low-, medium- and high-density regions). The meanings of the symbols are the same as in the left panel. 



bluer galaxies tend to have stronger Ha emission in an average 
sense. There are also some Ha emitters on the red sequence. These 
populations are very rare, but interestingly, we find that such Ha 
emitters with red colours are preferentially found in the cluster out- 
skirts (see below). On the other hand, MIR galaxies (shown as filled 
circles in Fig. UOt are more commonly seen on the red sequence. 
This is clearly seen in the histograms shown in Fig.QO] The Ha 
emitters are bluer and the MIR galaxies are redder. This may be 
because the MIR galaxies are dusty and red dened in spite of their 
strong star formation (see also lKovama et al.ll2008h . 

We also make the same plot for each environment (i.e. for low- 
, medium- and high-density regions) in the right-panels of Fig. [10] 
Interestingly, the red sequence appears quite differently with en- 
vironment. Most of the red galaxies in the high-density and low- 
density environment show little on-going star-forming activity. In 
contrast, a large number of red galaxies in the medium-density 
environment (i.e. in the cluster outskirts, groups and filaments) 
are detected in Ha and/or MIR. We detect 6 red Ha emitters in 
total, of which 5 are found in the medium-dens ity environment. 
We have already reported in lKovama et alj d2008l) a high fraction 
of MIR-detected dusty red galaxies in the medium-density envi- 
ronment. The above result on the Ha emitters supports this idea. 
ISaintonge. Tran. & Holdenld2008h showed that the number of MIR 
galaxies on the red sequence increases towards distant clusters up 
to z ~ 0.8. Our current result suggests that the MIR view of the 
red sequence is also dependent on environment at z ~ 0.8. 

6.2 Colour-colour plot 

We show in Fig.QT]the colour-colour diagrams (R—z v.s. z'—J). 
Using the three bands that neatly straddle the rest-frame 4000A 
break, we can effectively distinguish dusty red galaxies from pas- 
sively evolving red galaxies. This method is similar to the one 



used to distinguish between dusty and passive galaxies for the 
sample of the extremely r ed objects (EROs) at z ~ 1.5 (e.g. 
IPozzetti & Mannuccill2000]) . 

In the top-left panel of Fig.[TT] we notice that the star-forming 
galaxies (i.e. Ha emitters and MIR galaxies) are distributed in an 
elongated region along the direction of the reddening vector (shown 
by the arrow). In contrast, non-star-forming galaxies (small dots) 
are concentrated at around R — z ~ 1.6 and z — J ~ 0.6. Thus, 
we can separate between passive galaxies and star-forming galax- 
ies on this diagram. Many of the MIR galaxies are red in a single 
colour in Fig.QJJ but with an appropriate colour combination here, 
the passive galaxies and the dusty star-forming galaxies occupy dif- 
ferent areas on the diagram and can be well separated. 

In Fig.QT] we also make the same plot for each environment. 
As is also seen in Fig. QJJJ Ha emitters are blue and MIR galaxies 
are redder. Most of the red star-forming galaxies are detected only 
in MIR and even not in Ha. Such red MIR galaxies are frequently 
seen in the medium-density environment. We also find that very 
blue galaxies are detected only in Ha but not in MIR (as seen in 
Fig. [10}. This may be because bluer galaxies tend to be lower-mass 
galaxies (see Fig. UOt and to have lower SFR (see Fig. |6), so that 
we could not detect many of these low-mass star-forming galaxies 
by our MIR observation due to the MIR flux limit. 

6.3 Colour-density plot 

As we showed in the previous sections, the nature of star-forming 
activity of galaxies is related to environment (i.e. galaxy density). 
In Fig.[T2] we plot R — J colours as a function of the local density. 
This plot would be a good summary of our important findings. 

Firstly, colours of galaxies start to change at the medium- 
density environment. Clearly, the number of blue galaxies (with 
R — J < 2.0) decreases in high-density environment (log E 5 > 
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Figure 11. The colour-colour diagram (R — z' vs. z' — J) for all the 
observed fields (top-left panel), the high-density region (top-right), the 
medium-density region (bottom-left), and the low-density region (bottom- 
right), respectively. The small open circles represent all the member galax- 
ies. The open squares and the filled circles indicate the Ha emitters and 
the MIR galaxies, respectively. The arrows indicate reddening vectors cor- 
responding to E(B - V)=0.2. 
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Figure 12. The colour-density plot for all the member galaxies (open cir- 
cles). The small and large open squares indicate the weak (EW<50A) and 
strong (EW>50A) Ho emitters, respectively. The filled circles show the 
MIR galaxies. The dashed lines at log E5 = 1.65 and 2.15 indicate the di- 
viding lines between the low-, medium- and high-density regions as defined 
in Koyama et al. (2008). The solid-line shows the locus of the fraction of red 
galaxies calculated using our NIR-detected members. The vertical and hor- 
izontal error-bars indicate Poisson errors and the size of the environmental 
bins used for f(red), respectively. 



2.15). This is quantified by calculating the fraction of red galaxies 
(with R — J > 2.0) as a function of the local density (the solid-line 
locus in Fig. 1 12b. This res ult is consistent with our previous find- 
ing in lKovama et alj d2008l) . who made a similar plot and examined 
R — z' colour using all optical sources. We confirm the trend that 
galaxy properties are changed not only in cluster environment but 
also in relatively low-density groups and filaments. 

Secondly, we notice that most of the Ha emitters are blue, but 
there exist some red emitters with R — J > 2.0. These galaxies are 
candidates of dusty star-forming galaxies, and interestingly, these 
galaxies are preferentially found in the medium-density environ- 
ment. In fact, we have 6 such red Ha emitters in total and 5 out 
of 6 are in the medium-density environment. Also, 4 out of 6 such 
emitters are detected in MIR. This indicates that these red Ha emit- 
ters are not just ceasing their star-formation. Rather, it seems that 
they have very strong star-formation but are heavily attenuated by 
dust. 

Thirdly, the MIR galaxies are generally redder than the Ha 
emitters, and some of them are located on the red sequence. Such 
optically red MIR galaxies are also preferentially found in the 
medium-density region (i.e. transition environment) and some of 
them are not detected in Ha. 

It turns out that a fraction of the red galaxies in the cluster 
surrounding regions are dusty star-forming galaxies, which would 
have been difficult to be distinguished from passively evolving red 
galaxies with optical colour information alone. In contrast, most of 
the red galaxies in the high-density environment are passive galax- 
ies with little on-going star-forming activity. This may indicate that 
we are witnessing the transition phase of galaxies in the cluster out- 
skirts when/where some galaxies experience an active, but dusty 



star-forming phase before being truncated. More detailed discus- 
sion will follow in the next section. 



7 HIDDEN ACTIVITY IN DISTANT CLUSTERS 

7.1 Comparison between Ha and MIR activity 

As suggested in the previous section, a significant fraction of star 
formation activity is likely to be hidden in the optical data even 
in the Ha data, in particular in the surrounding regions of high- 
z clusters. To quantify the amount of hidden star formation, we 
calculate the ratio of SFR(IR) to SFR(Ha) for our MIR-detected 
galaxies down to 3.5a in detection (~ lOM /yr). Here, we mea- 
sure SFR(Ha) for the members that satisfy J-NB119>0.15 in- 
stead of J— NB119>0.3 that we used in Section 3.1, so that we can 
reach to fainter MIR-detected galaxies. In Fig. [T3] we plot the ra- 
tios, SFR(IR)/SFR(Ha), as a function of SFR(IR). Here, we do not 
apply for any extinction correction in deriving SFR(Ha). We see a 
weak trend that galaxies with higher SFR(IR) tend to have higher 
SFR(IR)/SFR(Ha) ratios. This is consistent with previous studies 
which show that the dusty component of star format ion rate be- 
comes higher as total star formation rate g oes up (e.g. Geac h et all 
l2006tlDopita et alj2002l ; lGarn et alj2009n . We can also notice that 
red MIR galaxies (shown as filled symbols in Fig. 113b tend to have 
higher SFR(IR)/SFR(Ha) ratios. This is reasonable because their 
red colours can be the results of heavy dust attenuation. 

The absolute values of SFR(IR)/SFR(Ha) span widely 
from ~2 to ~40, depending on their SFR(IR). For moder- 
ately star-forming galaxies (i.e. SFR(IR)~lOM0/yr), the ex- 
tinction of Ha strength seems to be ~ 1 mag, in excellent 
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agreement with the study of local star-forming galaxies (e.g. 
iKennicutt. Tamblvn. & Congdoii 1994). However, for actively star- 
forming galaxies (i.e. SFR(IR)> 2OM0/yr), the magnitude of ob- 
scuration is significantly large and An a can exceed ~3 mag in 
the extreme cases, even though Ha line is considered as one of 
the most reliable SFR indicators. There have been some pieces 
of evidence for such large amount o f obscuration in actively star- 
forming IR galaxies. For example, Poggianti & Wu (200Cj) pre- 
sented spectral natures of IR galaxies at intermediate redshifts, and 
showed that these galaxies tend to be classified as 'e(a)', which 
show emission lines as a sign of moderate star formation activity 
and strong balmer absorption at the same time. In fact, they showed 
SFR(IR)/SFR(Ha)~ 10-70, an d this is in good agr eement with our 
values. In cluster environments. ISmail et al] ( 1999) conducted deep 
radio observations in the CL0939 cluster at z = 0.41. They de- 
tected radio continuum light at 1.4GHz, as a sign of on-going star 
formation activity, from the galaxies which have been classified as 
post-starburst galaxies by optical spectra based on the presence of 
strong ba lmer absorptio n lines and the absence of emission lines. 
Recently, iDressler et alj d2009h conducted MIR observation of the 
A851 cluster at z — 0.41. They detected MIR emission not only 
from the majority of e(a) galaxies but also interestingly from ~30% 
of k(a) galaxies (i.e. post-starburst galaxies with no emission lines). 
These results all indicate that the optical signs of star-forming ac- 
tivity tend to be weak especially for galaxies with strong star for- 
mation that are well traced by MIR observations, and thus can be 
sometimes misleading. 

In fact, our MIR galaxies in RXJ1716 do not show strong 
Ha emissions. As shown in Fig. [6] we find that EW(Hq+[NII]) 
and SSFR(Ha) of the MIR-detected galaxies tend to be lower than 
those of the MIR-undetected Ha emitters. Also, the MIR galaxies 
are sometimes 'not' detected as Ha emitters. Unfortunately, we do 
not have any spectroscopic information for these galaxies yet. It 
is possible that these galaxies are even more heavily obscured ob- 
jects, but in the current situation, we cannot exclude the possibility 
that these galaxies are located at slightly lower or higher redshifts 
and their Ha lines do not fall within the NB119 filter. Therefore, 
we limit the sample to the galaxies with Ha detection in this paper 
because these galaxies are reliable members. We will present spec- 
troscopic data of this cluster to reveal the nature of MIR galaxies 
undetected in Ha in a forth-coming paper. 
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Figure 13. The ratios of SFR(IR) to SFR(Ho) plotted against SFR(IR). We 
plot only the Ha emitters with J colr — NB119X). 15 that are detected as 
single objects at MIR with detection significance of >3.5<r. The squares, 
triangles, and the circles indicate the galaxies in the low-, medium-, and 
high-density environments, respectively. The filled symbols show the red 
galaxies with R — J > 2.0. A representative error bar is shown as the 
cross. The uncertainty in SFR(Ho) corresponds to a Icr error in the aperture 
photometry, while we adopt a 50% error in SFR(IR) to take into account the 
large uncertainty in the conversion from vL v ^iim to -L(IR) (see text). 



and low-density environments, but no such galaxies are found in 
the high-density region. These results are consistent with our find- 
ings in Section 6 where we find many dusty red galaxies in the 
medium-density regions. Such enhanced dusty star-formation in 
the intermediate-density environme nt has also been reported for 
lower- z samples in the literature (e.g. Wolf, Gray, & Meisen heimeJ 



l2005tlFadda etalj200 8: Gallazz i et al.l2009l ; iTran et al . 2009). Out- 
current results show a similar trend up to z ~ 0.8. Our unique 
wide-field Ha and MIR surveys in and around the RXJ1716 cluster 
demonstrate the necessity of wide-field IR observations covering 
the transition environment in the cluster suburbs in order to reveal 
the key processes to shape the environmental dependence of galaxy 
evolution and its link to IR activity of galaxies. 



7.2 Environmental dependence of the hidden activity ? 

Taking the unique advantage of the wide-field coverages of our 
both Ha and MIR imaging surveys presented in this paper, we 
attempt to examine the environmental dependence of 'hidden' 
star-formation around a z ~ 0.8 cluster for the first time. In 
Fig. [13] we use different symbols depending on their environment 
(squares for low-density, triangles for medium-density and circles 
for high-density). Interestingly most of the 'heavily attenuated' 
galaxies (e.g. SFR(IR)/SFR(Ha)^15) are located in the medium- 
density regions. Although the number of MIR galaxies is not large, 
there is a hint that very dusty galaxies most frequently appear in 
groups/filaments where galaxy properties are dramatically changed 
(see Section 6.3). Note that some of those galaxies in the outskirts 
may have large relative velocities with respect to the cluster core 
and their Ha lines may fall on the wings of the transmission curve 
of the NB119 filter, and therefore their SFR(Ha) may be under- 
estimated. Their red colours, however, indicate that they are in- 
deed dusty objects. Also, it can be seen that bright MIR galax- 
ies (e.g. SFR(IR) S; 40Mo/yr) are found only in the medium- 



8 GLOBAL EVOLUTION OF STAR-FORMATION 
ACTIVITY IN CLUSTERS 

8.1 Ha view of the cluster evolution 

In this final section, we discuss the global evolution of star for- 
mation activity in clusters together with the previous Ha cluster 
studies in the literature. iFinn et all d2005h sketched the evolution of 
Ha-derived total SF in clusters (SSFRHa) and that divided by clus- 
ter mass (SSFRH a / M e i) as a function of cluster redshift (see also 
Kod ama et ai] |2004). We now estimate the se quantities (SS FRhq 
and £SFR Ha /M ci ) for the RXJ1716 cluster. lFinn et al.1 J2005h used 
all the cluster member galaxies within 0.5 x -R200 for each cluster 
to derive ESFRH a and divide ESFRhq by M200 for each cluster 
to derive ESFRhh/A/c!. The radius and mass (-R200 and M200) of 
each cluster are defined by the following equations: 

^- 2 - 47 iouoL/ s V n A + » o( i + ,)3 Mpc (4) 
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1.71 x 10 1 



Vioookm/ s ; ^/n A + n (i + z ) 



-.M a 



(5) 

where a is the velocity dispersion of the cluster. However, for 
an unrelaxed cluster such as this RXJ1716 cluster, the veloc- 
ity dispersion tends to be over-estimated due to the presence of 
kinematical substructures. In fact, the velocity disper sion of the 
RXJ17 16 cluster is estimated to be ~1500km/s in iGioia et all 
(1999), which is much larger than that estimated from its X- 
ray temperature (~6 00km/s) as t h ey no ted. Based on the new X- 
ray measurement in lEttori et al.l ( 120040 . we re-estimate the ve- 
locity dispersion to be a =837 km/s. Based on this value, we 
derive #200=1-3 Mpc and M 20 o=6.5 x 10 14 M© for the RXJ1716 
cluster. We then integrate SFR(Ha) of the member galaxies de- 
tected at more than 3cr level in Ha within 0.5xi?200, and correct 
for incompleteness due to the wings of the filter response func- 
tion by multiplying a factor ~1.4 to derive ESFRh q . We show 
our measurements, ESFR Ha = 210 ± 42 Mp /yr, in Fig. [0] 
with t hose of other clus ters in iFinn et al] fc005h (i.e. orig inally 
from IBalogh et alj |2002| fo r A1689, IBalogh & Morris! l200d for 
A2390, ICouch et al] l200l1 for AC 114. iKodama et al.l [2004, for 
l2005l for CL1040, CL1054-12 and CL1216, 



CL0024, 



Finn ct al 



and IFinn. Zaritskv. & McCarthvl |2004| for CLJ0023+0423B). As 
can be seen in Fig. [14] we find no significant evolution in £SFRh q 
with redshift (lower panel), while £SFRh q /M c ( seems to increase 
towards higher redshift clusters (upper panel). Although the scatter 
at a given redshift is very large, ESFRhq/A/c; approximately scales 
with (1 + z) 6 as shown by a solid curve in the top panel of Fig. 1141 
It is possible that the scatter among clusters at a similar 
redshift is due to the difference in cluster mass. As shown in 
Fig. E] ES FRrfa/M e i i s wel l correlated with cluster mass, in con- 
sistent with IFinn et alj d2005h (see also IPoggianti et alJPioO^ and 
iHomeier et alj|2005h . It may indicate that the cluster mass is also 
an important factor to determine cluster properties. We should note, 
however, that the determination of cluster mass is highly uncer- 
tain and this could be the biggest source of uncertainty in drawing 
any general picture of cluster evolution. In fact, if w e take a ~ 
1500k m/s as derived from the optical spectroscopy in IGioia et alj 
dl999r> . ESFRrfa increases by a factor of ~2 and ESFR Ha /M ci de- 
creases by a factor of ~ 3. We would need a much larger sample of 
clusters to obtain more secure trend. 



8.2 MIR view of the cluster evolution 

In the previous sub-section, we have examined the evolution of star- 
forming activity in clusters based on Ha intensity assuming a con- 
stant 1 -magnitude extinction correction for all the emitters. This 
assumption is valid for moderately or weakly star-forming galax- 
ies, but much larger extinction correction would be needed for very 
actively star-forming galaxies (see Fig. 1131). In order to take into 
account the hidden star-formation as much as possible, we now 
use SFR(IR) for the MIR-detected galaxies while we stick to use 
SFR(Ha) with 1-mag correction for the MIR-undetected galaxies, 
and sum up all the SFR within 0.5 x #200 from the cluster centre. 
In this calculation, we included the blended MIR sources as well, 
and used the combined MIR flux of such sources. These all sum 
up to the total SFR(Ha+MIR) of ~ 620 ± 250 M Q /yr , which is 
more than factor 2 larger than the integrated SFR(Ha) even after 
applying for the constant 1-mag extinction correction. Note that we 
adopt a 50% uncertainty in ESFRir, which is the same amount as 
what we assumed in SFR(IR) for individual galaxies (Section 2.3). 
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Figure 14. (bottom): The integrated total SFR(Ho) within 0.5xR 2 oo as a 
function of redshift. (top): The integrated total SFR(Hq) normalized by the 
cluster mass Af2oo as a function of redshift. 
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This may be an overestimation for an integrated quantity ESFRir, 
but we take it as a conservative error. 

We suggest that non-negligible fraction of star-forming ac- 
tivity in clusters is hidden by dust and thus strongly supports the 
importance of IR study of distant clusters and/or accurate extinc- 
tion correction to obtain true SFRs and to quantify the evolution 
of star-forming activity in clusters. Although the number of known 
clusters studied so far in MIR is very small, some authors have 
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attempted to sketch the global evolution of st ar formation activ- 
ity in clusters b ased o n the MIR-derived S FRs ( Geach et al l 20061 : 
iBai et alJl2007L 2009; lHaines et alJl2009al: iKrick et al.ll2009h . We 
here follow the scheme in IBai et alj faOVftL who measured inte- 
grated cluster SFR(IR) within 0.5xi?200 from the cluster centre. 
Their flux limit in IR corresponds to ~lOM0/yr at z ~ 0.8, but 
they integrated SFR down to ~2A/©/yr by extrapolating their IR 
luminosity function to that depth. We can go down to the sim- 
ilar depth by combining the MIR- and Ha-derived SFRs within 
0.5x7?200 of the RXJ1716 cluster. In Fig. [T6l we show the total 
SFR(Ha+MIR) normalized by clu ster mass for R XJ1716 cluster 
with SFRs(MIR) of ot her clusters inlBai et al.1 120071) (i.e. originally 
from MIR surveys by Bai et al 120061 f or Coma, Biviano et al]2004l 
forA2218 l Ducetalj2002l for A1689.ICoiaetal.l2005al forA2219. 
ICoia et al.ll2005bl f or C10024, iMarcillac et al]|2007l for RXJ0152 
and lBai et alj|2007l for MS 1054). We also add th ree more recentl y 
reported data points at z = 0.06 (A3266) fro m lBai et alj j2009t> . 
at z = 0.08 (A2255) from lShim et~ai] 120091) and at z = 1 as a 
compo site value of three z ~ 1 clusters presented in lKrick et al.1 
(2009). We use M200 as the cluster mass for RXJ1716 (same as in 
Section 8.1), but the trend does not change if we u se a lensing mass 
of 2. 6 ± 0.9/1" 1 x 1 14 M Q Iciowe et alJl 19981) to be consistent 
with lBai et alj 120071) who mainly use lensing masses. We can see 
a similar evolutionary trend with redshift following oc (1 + z) 6 in 
cluster star-forming activity to that seen in Fig. 1141 

On the other hand, many field studies showed that the cosmic 
star formation rate density or SSFR of field galaxies decline since 
z ~ 1 following tx (1 + z) 3 (e.g. lHoDkinsll2004 lYoshida et al.1 
2006). The difference in the power index between cluster and field 
suggests that the galaxy evolution in cluster environment is accel- 
erated compared to that in the general field since z ~ 1 (see also 
Kodama & Bower 2001), although we need a larger sample of dis- 
tant clusters to confirm this interesting result and to quantify envi- 
ronmental effects. 

It would be expected that the number (or fraction) of actively 
star-forming galaxies and the relative importance of hidden star- 
formation may increase with redshift. The combined approach of 
deep Ha and MIR observations for z » 1 clusters is thus essen- 
tial and it will enable us to quantify the importance of hidden side 
of star formation activity in clusters and its evolution. 



9 SUMMARY 

We conducted a narrow-band Ha imaging survey for the RXJ1716 
cluster at z — 0.81 with MOIRCS on the Subaru Telescope. This is 
the first wide-field Ha emitter survey for z > 0.8 clusters to date, 
and we neatly covered the pre-defined filamentary large-scale struc- 
tures in and around this cluster. Combining with the wide-field MIR 
imaging data taken with AKARI satellite, we mapped out not only 
normal star-forming galaxies but also dusty starbursting galaxies 
across a wide range in environment. Our findings are summarized 
as follows: 

(i) The spatial distribution of the Ha emitters are very similar 
to that of the MIR galaxies. We find that the Ha emitters avoid 
cluster central region (R c <0.25Mpc) just as the MIR galaxies do, 
and that the fraction of Ha emitters is higher in the lower-density 
environment. 

(ii) Ha emission lines are detected mainly from blue galax- 
ies, while the MIR galaxies tend to be redder. There are some 
MIR galaxies on the red sequence, and interestingly, such red 
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Figure 16. The integrated total cluster SFR(IR) per unit cluster mass as 
a function of redshift. The open circles are taken from Bai et al. (2007). 
The open square and the diamond are from Bai et al. (2009) and Shim et 
al. (2009), respectively. The open triangle shows the composite value of 3 
clusters at z = 1 in Krick et al. (2009). For RXJ1716, we combine the Ha 
data and the MIR data to derive the total SFR (see text). 



MIR galaxies are most commonly seen in the 'medium-density' 
environment such as cluster outskirts, groups and filaments. We 
also find that the MIR galaxies are not strong Ha emitters. Also, 
there are some MIR galaxies without detectable Ha emission lines 
(EW^30A). These suggest the existence of large amount of highly 
obscured galaxies in the distant cluster environment. However, we 
still need spectroscopic confirmation of membership of the individ- 
ual MIR and Ha sources to prove this interesting result. 

(iii) We find some Ha emitters on the red sequence. Although 
such population is rare (only 6 in total), these red Ha emitters 
are concentrated in the medium-density environment, similar to the 
spatial distribution of MIR sources. We also find that many of the 
red Ha emitters (4 out of 6) are detected in MIR. This suggests that 
such red emitters are not just gradually truncating their star-forming 
activity but actively star-forming with heavy dust obscuration. 

(iv) We examined the amount of hidden star formation based on 
the ratio of SFR(IR)/SFR(Ha) for the MIR galaxies. We find that 
for moderately star-forming galaxies with SFR(IR)< 10Mo/yr the 
amount of extinction for Ha is ~1 mag, in excellent agreement 
with local spirals. However, for actively star-forming galaxies with 
SFR(IR)>2OAf0/yr, Ha lines are more heavily attenuated and the 
extinction can exceed Ah q ~3 in the extreme cases. Interestingly, 
most of such very dusty galaxies are located in the medium-density 
environment, suggesting the enhancement of star formation activity 
in such environment, most likely triggered by galaxy-galaxy inter- 
actions and/or mergers. 

(v) Combining our unique Ha and MIR data, we derive total 
SFR in the cluster down to ~2A/o/yr within 0.5xi?200 from the 
cluster centre. We find that the cluster total SFR based on Ha alone 
can be underestimated about factor ^2 even after 1-mag extinc- 
tion correction. We suggest that the mass-normalized cluster SFR 
evolves with redshift following tx (1 + z) 6 , although the number 
of clusters studied so far is still small. We need a much larger sam- 
ple of clusters studied in MIR and should go back to more distant 
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clusters at the same time, to fully understand the evolution of star- 
forming activity in clusters. 

Overall, we have stressed the importance of cluster outskirts 
as a key environment for galaxy transition. Our results also imply 
that we need to be very careful when studying star formation ac- 
tivity in distant clusters. We may miss a significant amount of star 
formation if we have only optical information (even with the Ha 
line), in particular in the cluster outskirts where starbursting popu- 
lations are seen. In order to trace the true star formation activity in 
distant clusters, MIR-FIR observations will be essential on top of 
the optical-NIR observations. 
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